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The thermal conductivity n of the iron-arsenide superconductor Bai_ a; K :! ;Fe2As2 (T c ~ 30 K) was 
measured in single crystals at temperatures down to T ~ 50 niK (~ T c /600) and in magnetic fields 
up to H — 15 T (~ _ff C 2/4). A negligible residual linear term in k/T as T — > shows that there 
are no zero-energy quasiparticles in the superconducting state. This rules out the existence of line 
and in-plane point nodes in the superconducting gap, imposing strong constraints on the symmetry 
of the order parameter. It excludes d-wave symmetry, drawing a clear distinction between these 
superconductors and the high-T c cuprates. However, the fact that a magnetic field much smaller 
than H C 2 can induce a residual linear term indicates that the gap must be very small on part of the 
Fermi surface, whether from strong anisotropy or band dependence, or both. 

PACS numbers: 74.25.Fy, 74.20.Rp,74.70.Dd 



Several experiments have been performed to address 
the symmetry of the order parameter in iron arsenide su- 
perconductors. Early studies on mostly polycrystallinc 
samples of i?FeAs(0,F), where R = La, Nd or Sm, 
have led to indications that appear contradictory. While 
point-contact spectroscopy [ljj, ARPES Q and penetra- 
tion depth [3], 0, 0] point to a full superconducting gap 
without nodes, specific heat [|| and NMR 0] data were 
interpreted in terms of a nodal superconducting gap. Re- 
ports on single crystals of doped BaFe2As2, also appear 
contradictory. In (Ba,K)Fe2As2, ARPES studies have 
found an isotropic superconducting gap with a magni- 
tude of 12 meV on one Fermi surface and 6 meV on an- 
other 0, 0, [l(| , specific heat measurements are broadly 
consistent with an s-wave gap of 6 meV 111! , as are muon 
measurements of the supcrfluid density [12J . By contrast, 
penetration depth studies find a power-law variation [l3| , 
as in Co-doped BaFe2As2 [3, as opposed to the expo- 
nential temperature dependence expected of an isotropic 
s-wave gap. 

In an attempt to shed further light on the structure 
of the superconducting gap, we have measured the ther- 
mal conductivity k of K-doped BaFe2As2. Heat trans- 
port is a powerful probe of symmetry-imposed nodes in 
the superconducting gap [l5| . We find a negligible resid- 
ual linear term kq/T in k/T as T — > 0, strong evidence 
that there are no nodes in the gap of this superconduc- 
tor. Indeed, a line of nodes would have given a sizable 
and universal (i.e. impurity independent) kq/T flBl ], as 
in cuprates (lfjj , ruthenates (l7j and some heavy- fermion 
superconductors [l8| . Given the large impurity scattering 
rate in our samples, point nodes would also have given 



a sizable (albeit non- universal) kq/T [15[ - unless they 
happend to lie along the c-axis, perpendicular to the di- 
rection of heat flow in our measurements. However, a 
magnetic field H applied along the c-axis induces a finite 
kq/T even when H << H C 2- This shows that the super- 
conducting gap must be very small on some part of the 
Fermi surface, either because of a pronounced anisotropy 
on one Fermi surface (whereby the gap has a deep min- 
imum in some directon) or because of pronounced band 
dependence, causing one Fermi surface sheet to have a 
very small gap. 

Single crystals of Bai_ K K l: Fe2As2 were grown from 
FeAs flux [19( . The doping level of the two samples used 
in the present study, labeled A and B, was determined 
from the c-axis lattice parameter, giving x= 0.25 for sam- 
ple A (T c = 26 K) and x= 0.28 for sample B (T c = 30 K), 
both on the underdoped side, i.e. below optimal doping 
(x ~ 0.4) [i^ . 

Samples were cleaved into rectangular bars with typi- 
cal size 1.5 x 0.3 x 0.05 mm 3 . Silver wires were attached 
to the samples with a silver-based alloy, providing ultra- 
low contact resistance of the order of 100 Thermal 
conductivity was measured along the [100] direction in 
the tetragonal crystallographic plane in a standard one- 
heater-two thermometer technique (2lj . The magnetic 
field H was applied along the [001] tetragonal axis. All 
measurements were done on warming after cooling in con- 
stant H from above T c to ensure a homogeneous field 
distribution in the sample. All aspects of the charge and 
heat transport are qualitatively the same in both sam- 
ples, with minor quantitative differences. For simplicity, 
only the data for sample A is displayed here. Whenever 
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FIG. 1: (Color online) Temperature dependence of the in- 
plane resistivity p(T) of Bai_ a; K; r Fe2As2, with x ~ 0.25 and 
T c = 26 K. Inset: same data plotted as a function of T 2 . The 
line is a linear fit to p(T) = po + AT 2 below 70 K. 



quoted, we provide quantitative values for both samples. 

In Fig.[Tl we show the electrical resistivity p(T) of sam- 
ple A as a function of temperature. Below ^150 K, 
p(T) shows a notable downturn, followed by a range 
where p(T) is well described cither by a power-law de- 
pendence, p(T) = po + AT n , with n between 1.6 to 1.8, 
or by p(T) = p Q + AT 2 , below 70 K or so. Using the 
latter fit, we get po — 47 pil cm (sample A) and 28 
pfl cm (sample B). Since the same contacts are used for 
electrical and thermal transport, this allows us to ac- 
curately estimate the normal-state thermal conductivity 
k-^/T in the T — > limit, via the Wiedemann- Franz law, 

kn/T = Lo/po where L = ^-( ! f) 2 , giving k n /T = 520 
pW/K 2 cm (sample A) and 875 pW /K 2 cm (sample B). 

Residual linear term in zero magnetic field. The ther- 
mal conductivity k(T) of sample A, measured for H=0, is 
shown in Fig. El The data is plotted as k/T vs T 165 . The 
linear fit displayed in Fig. [2] shows that the data below 
0.4 K is well described by the function k/T = a + bT a , 
with a = kq/T = 5 pW / K 2 cm and a = 1.65. This same 
function describes the data of sample B equally well, with 
a = K a /T = 7 pW I K 2 cm and a = 1.5. The first term is 
the residual linear term of interest here [IH . The second 
term is due to phonons, which at low temperatures are 
scattered by the sample boundaries. Although the latter 
term would be expected to give k p oc T 3 , i.e. a = 2, mea- 
surements on single crystals with smooth surfaces rarely 
show this textbook behavior because of a specular reflec- 
tion off the surface, and in practice 1 < a < 2 mm hi. 

The magnitude of the residual linear term extracted 
from the fits in Figgis extremely small. In similar mea- 
surements on samples where no residual linear term is ex- 
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FIG. 2: (Color online) Temperature dependence of the ther- 
mal conductivity k(T), plotted as k/T vs T 1 ' 65 , in zero mag- 
netic field. The line is a linear fit over the temperature range 
shown. 



pected, Ko/T was indeed found to be zero, within an error 
bar of ±5 pW / K 2 cm or so 

mill- Within those error 
bars, our two samples of Bai_ 2; K a ,Fe2As2 exhibit negligi- 
ble residual linear terms. Let us put these minute kq/T 
values into perspective. Comparison with the normal- 
state conductivity, gives the ratio (kq/T)/(k^/T) ~ 1% 
in both samples. Secondly, these kq/T values are much 
smaller than theoretical expectations for a nodal super- 
conductor. (For a gap without nodes, kq/T should be 
zero [HI].) For a quasi-2D d-w&ve superconductor, with 
four line nodes along the c-axis, the residual linear term is 
given, in the clean limit (hT « A ), by [HiESI^I 
kq/T = (fcj|/6c)(fc,FUF/Ao), where c is the interlayer sep- 
aration, kp and vf the Fermi wavevector and velocity at 
the node, respectively, and Ao the gap maximum. Tak- 
ing c = 6.6 1, v F = 0.3 eV 1 = 0.45 x 10 5 m/s and 
a Fermi wavevector kp = 0.35 xn/a = 0.28 A -1 , and as- 
suming a weak-coupling Ao = 2.14 k^T c , we get kq/T = 
140 pW I K 2 cm. This is at least 20 times larger than the 
value extracted from our fits to the data. In those ma- 
terials where universal heat transport has been verified, 
the measured value of kq /T is in good quantitative agree- 
ment with this theoretical expectation 0,111, Thus 
we can safely conclude that the gap in Bai_ 2 ,K 2; Fe2As2 
does not contain a line of nodes anywhere on the Fermi 
surface. In particular, this rules out d-wavc symmetry, 
whether d x i_ y i or d xy . This result sets cuprates and iron 
arsenides apart as two distinct types of high-temperature 
superconductors. 

It is in principle possible for the superconducting gap 
in pnictides to have point nodes as opposed to line nodes, 
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FIG. 3: (Color online) Left panel: Temperature dependence 
of the thermal conductivity measured in magnetic fields of 0, 
1,4, 10 and 15 T (from bottom to top), plotted as k/T vs T 2 . 
Solid lines are a linear fit to each curve in the range shown. 
Vertical dashed lines indicate T = 0.1 K (blue) and T = 
0.15 K (green). The value of k/T at those two temperatures 
is plotted vs magnetic field in the right panel. Right panel: 
isotherms of k/T as a function of magnetic field H , for T — > 
(obtained by extrapolating the linear fits in the left panel 
to T = 0), 0.1 K and 0.15 K. In all three cases, k/T rises 
approximately linearly with H, with the same slope. The 
solid line is a linear fit to the T — > data, also reproduced in 
Fig.H 



one of the scenarios suggested by penetration depth stud- 
ies [3]. The zero-energy quasiparticles associated with 
point nodes give a residual linear term which grows with 
impurity scattering (26[, so that kq/T can become a 
substantial fraction of k^/T (26[. Given po, we can 
estimate the normal-state impurity scattering rate Tq 
roughly from the plasma frequency uj p = c/Ao, where 
Ao is the penetration depth, approximately equal to 200 
nm 0. This gives T /T c ~ 1.7 and 0.9 for sample A 
and B, respectively. This is very substantial, and would 
give a large residual linear term, comparable to the case 
of the line node. We therefore conclude that point nodes 
in the gap of Bai_ K K l: Fe2As2 are also unlikely, unless 
they are located along the c-axis and do not contribute 
to in-plane transport. Needless to say, our data also rule 
out the possibility of an entirely gapless (or ungapped) 
Fermi surface, proposed by some authors [2^], at least 
down to the 1 % level. 

The nodes we have discussed so far are imposed by 
symmetry, the result of a sign change in the order param- 
eter around the Fermi surface. Such symmetry-related 
nodes are broadened by impurity scattering, giving rise 
to a sizable Ko/T. The superconducting gap can also 
go to zero in certain directions because of a pronounced 
anisotropy that is not imposed by symmetry. However, 
such "accidental" nodes in a gap with s-wave symme- 
try will be lifted by impurity scattering, making the gap 
more isotropic [3(| . Our zero-field data is consistent with 
an s-wave gap, including one with strong anisotropy. 

Field dependence of thermal conductivity. The effect of 
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FIG. 4: (Color online) Residual linear term kq/T in the ther- 
mal conductivity of various superconductors as a function of 
magnetic field H, presented on normalized scales k 3 /kn vs 
H/H C 2- For Bai_ I K 2 ;Fe2As2 (red circles), we use km /T ob- 
tained from the Wiedemann-Franz law applied to the extrapo- 
lated residual resistivity po (see text) and H C 2=65 T [3l| . The 
data for the clean and dirty isotropic s-wave superconductors 
Nb and InBi, respectively, are reproduced from [24|]. The 
data for the d-wave superconductor Tl-2201 is for a strongly 
overdoped sample (T c = 15 K and H c2 ~ 7 T) The 
data for multi-band superconductor NbSe2 (T c = 7 K and 
H c2 = 4.5 T) is from [2j|. 



a magnetic field H (H || c) on the thermal conductivity 
of Bai_ 2; K a ;Fe2As2 is displayed in Fig. [3l In the panel 
on the left, k/T curves are seen to shift upwards almost 
rigidly with field. In the right panel, the value of k/T at 
three temperatures (T -> 0, T = 0.1 K and T = 0.15 K) 
is seen to rise linearly with H up to our highest field of 
15 T, which corresponds roughly to -ff C 2/4 (using H c2 at 



T — » as approximately 65 T [3l|). In Fig. [4] we com- 
pare the field dependence of Ko/T in Bai_ a; K 2 ;Fe2As2 
with the dependence in various other superconductors, 
using normalized conductivity and field scales, k s /kn 
and H/H C 2- In a d-wave superconductor like the over- 
doped cuprate Tl-2201 (T c = 15 K and Hr3~ 7 T), k /T 
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rises very steeply at the lowest fields 
as kq/T oc s/H, following the density of delocalized zero- 
energy quasiparticles outside the vortex cores j^. By 
contrast, in an isotropic s-wave superconductor like Nb, 
the rise in kq/T is exponentially slow at low fields as it 
relies on the tunneling of quasiparticles between localized 
states inside adjacent vortex cores, which at low fields are 
far apart. In many materials, however, the situation is 
not so clearcut. A good example are the multi-band su- 
perconductors MgB2 34 1 and NbSe2 [25| . Here the mag- 
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nitude of the s-wave superconducting gap is very different 
on two sheets of the Fermi surface. In both materials, the 
small gap is roughly one third of the large gap, so that a 
field H ~ H c2 /9 is sufficient to kill superconductivity on 
the small-gap Fermi surface, which can then contribute 
its full normal-state conductivity even deep inside the 
vortex state. Specifically, at H = H C 2/5, kq/T is already 
half (one third) of n N /T in MgB 2 (NbSe 2 ). It is one 
tenth in Bai_ a; K a; Fe2As2. By comparison, kq/T is still 
negligible in a single-gap superconductor like pure Nb or 
disordered InBi (see Fig. [J]). This shows that the super- 
conducting gap must be small on some part of the Fermi 
surface of Bai_ a; K a ;Fc2As2; relative to the gap maximum 
which controls H C 2- 

In Fig. IU we show data for NbSe2 [HI, where we see 
that Ko/T = a,t H = and kq/T rises linearly above 
H C 2/30, with a slope of 1.67 in the normalized units of 
Fig. |U In Bai_ z: K 2 ;Fe2As2, Ko/T also rises linearly, with 
a normalized slope of 0.46 in sample A and 0.21 in sample 
B. In a multi-band scenario, the magnitude of this slope 
is roughly proportional to the value of the normal-state 
conductivity kn/T of the small-gap Fermi surface rela- 
tive to the overall conductivity. The fact that the slope in 
absolute units is larger in sample A (3.7 /zW / K 2 cm T) 
than in sample B (2.8 iiW / K 2 cm T) eventhough its to- 
tal normal-state conductivity kn/T = Lo/po is smaller 
(520 vs 875 iiW/K 2 cm) is suggestive of a multi-band sit- 
uation with the impurity scattering rate being different 
on different Fermi surfaces. 

A recent heat transport study of the low-T^ nickel- 
arsenide superconductor BaNi2As2 (T c = 0.7 K) [35( gave 
K /T = 0, as here in Bai_ x K a; Fe2As2 (T c = 26 - 30 K). 
However, it found a much slower increase of kq/T at low 
H/H c 2, consistent with a gap that is large everywhere on 
the Fermi surface. 

We conclude that there are no nodes in the super- 
conducting gap of Bai_ a; K a ;Fe 2 As2, at least at x ~ 
0.25 — 0.28, with the possible exception of point nodes 
along the c-axis. This excludes (f-wave symmetry, and 
any other symmetry that requires line nodes on the multi- 
sheet Fermi surface of this superconductor. Symmetries 
consistent with this constraint include s-wave and s±, 
whereby a full gap changes sign from the electron Fermi 
surface to the hole Fermi surface (36| . From the rapid rise 
of kq/T with magnetic field at very low fields, we infer 
that the gap must be very small on some portion of the 
Fermi surface. This fc-dependence of the gap magnitude 
can come from angle dependence or band dependence, or 
both. In many experiments, the presence of a very small 
gap could mimic that of a node. 
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